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ABSTRACT 



Context. The Jacobi virial equation is a very powerful tool in exploring several aspects of the stellar internal structure and evolution. 
In a previous paper we have shown that the function [o/6] GR /A°- 9 (K) is constant (as 0.4) for pre main-sequence stars (PMS), white 
dwarfs (WD) and for some neutron star (NS) models, where aoR and Pgr are the form-factors of the gravitational potential energy 
and of the moment of inertia. 

Aims. To investigate the structural evolution of another type of celestial bodies, we extend these calculations to gaseous planets. We 
also analyse the cases for which this function is not conserved during some stellar evolutionary phases. Concerning NS, we study the 
influence of the equation of state (EOS) on this function and refine the exponent of the auxiliary function A(R). We also present a 
macroscopic criterion of stability for these stars. 

Methods. Non-stop calculations from the pre main-sequence to the white dwarf cooling sequences were performed with the MESA 
code. The covered mass range was 0.1-1.7 M Q . We used the same code to compute models for gaseous planets with masses between 
0. 1- 50 Mj. Neutron star models were computed using two codes. The first one is a modified version of the NSCool/TOV subroutines. 
The second code is a plain TOV solver that allows one to use seven previously described EOS. The relativistic moment of inertia and 
gravitational potential energy were computed through a fourth-order Runge-Kutta method. 

Results. By analysing the internal structure of gaseous planets we show that the function [ff/3] GR /A S (R) = T(M, EOS) is conserved 
for all models during the whole planetary evolution and is independent of the planet mass. For the pre main-sequence to the white 
dwarf cooling sequences, we have found a connection between the strong variations of T(M, EOS) during the intermediary evolution- 
ary phases and the specific nuclear power. A threshold for the specific nuclear power was found. Below this limit this function is 
invariant (w 0.4) for these models, i. e., at the initial and final stages (PMS and WD). For NS, we showed that the function T(M, EOS) 
is also invariant (» 0.4) and is independent of the EOS and of the stellar mass. Therefore, we confirm that regardless of the final 
products of the stellar evolution, NS or WD, they recover the initial value of T(M, EOS) ~ 0.4 acquired at the PMS. Finally, we have 
introduced a macroscopic stability criterion for neutron star models based on the properties of the relativistic product [a/J] GR . 

Key words, stars: evolution - stars: interiors - stars: pre-main sequence - stars: neutron - white dwarfs 
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£^ 1. Introduction constancy of [q^Jn^ by using realistic stellar models, Claret 

5sr | & Gimenez (1989) provided values of a and fi for several grids 

^ J The vectorial equations often used in physics are sometimes too of stellar evolutionary models. The availability of these form- 

. . complicated to be solved in a closed form. One of the advan- factors allowed deeper investigations into the Jacobi dynamics 

k*" fages of the virial theorem is that it reduces the complexity of of stellar evolution (Quiroga & Claret 1992ab; Quiroga & Mello 

°k> these equations, but this frequently comes at the expense of los- 1992, Quiroga & Cerqueira 1992). Unfortunately, due to limi- 

ing information. The balance is, however, positive because some tations in the opacity tables available at that time, these stellar 

r$ complicated problems can be treated and important information models only covered the pre main-sequence, the main sequence 

can be derived. The moment of inertia and the gravitational po- and only the first stages of the red giant branch. Therefore, many 

tential energy are fundamental elements that are handled by the important aspects of the later stellar evolution could not be inves- 

Jacobi virial equation. This equation is a very useful tool for ex- tigated in the mentioned papers. Another considerable limitation 

ploring several aspects of the stellar internal structure and evo- of these papers concerns mass loss, which is not considered. As 

lution. Given the nature of this equation, a mathematical rela- we will see in the next sections, this point is very important, 

tion between the moment of inertia and the gravitational poten- For example, a model with an initial mass of 1 .7 M Q at the pre 

tial energy reduces its complexity, which facilitates managing main-sequence reaches the white dwarf stage with only 0.6 M . 

it. Ferronsky et al. (1978) argued that this relationship could Modern stellar models available in the literature improved the 

be expressed by the constancy of the product a/3 in the Newto- situation: Claret (2006, 2007) published grids covering a wide 

nian approximation, where a and /3 are the form-factors of the range of chemical compositions (Z= 0.001-0.10), masses (0.8- 

gravitational potential energy and of the moment of inertia, re- 125 M Q ) and also offered detailed tracking of the late phases of 

spectively. However, the laws of mass distribution used by these stellar evolution. 

authors to evaluate the product [Qj6] Newt covered only a small In a recent paper Claret (20 12) showed that the product of the 

branch of the Herzsprung-Russel (HR) diagram. To check the form-factors a and /3 - neither the Newtonian nor the relativistic 
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one - is not conserved during the core hydrogen-burning phase, 
helium-burning, thermally pulsating asymptotic giant branch, or 
'blue loops' for models evolving from the pre main-sequence to 
the white dwarfs stages. In the same paper, it was shown that 
for very compact objects (neutron stars (NS) - for example) the 
effects of general relativity must be taken into account and the 
product a/3 computed using the Newtonian approximation must 
be superseded by the relativistic product [o'/SJqr. 

In our formulation the conserved quantity at the begin- 
ning and at the end of the stellar evolution is the function 
[a/3] CR /A 0S (R) = T(M,EOS) * 0.4. The auxiliary function 

A(R) is given by [l - 2G ^ 2 (g) J and is evaluated at the surface 
of each model, c is the speed of light, R is the radius, M(R) 
is the gravitational mass, and G is the constant of gravitation. 
Claret (2012) has also shown that, regardless of the final prod- 
ucts of stellar evolution (white dwarfs (WD) or NS), they recov- 
ered the fossil value of [oy6] GR /A 09 (R) acquired at the pre-main 
sequence (PMS) phase. The refinement of the exponent in the 
quoted function using an extended set of more elaborate EOS is 
given in Sect. 5. 

In the present paper we explore some aspects connected to 
these investigations. First, we extend the calculations of the mo- 
ment of inertia and the gravitational potential energy to gaseous 
planets and show that the function [a/3] CR / A & {R) (hereafter 
F(M, EOS)) is conserved during the planetary evolution, regard- 
less of the initial mass. Second, we investigate why this func- 
tion is not conserved during some evolutionary phases when we 
consider a complete evolution from the PMS to the WD cool- 
ing sequences (hereafter PMS-WD models). Third, we confirm 
that F(M, EOS) is conserved for NS and, in addition, is indepen- 
dent of the equation of state and of the stellar mass. We also 
derive, heuristically, an alternative criterion of stability based on 
the properties of the gravitational potential energy and of the 
moment of inertia. 



2. Neutron star, pre main-sequence, white dwarf, 
and gaseous planets models 

2.1. Planets 

The models of gaseous planets were computed using the MESA 
code (Paxton et al. 201 1, version 4298) for a chemical composi- 
tion of X=0.73 and Z=0.01. The adopted mixing-length param- 
eter omit is 1.5. We computed planet models from 0. 1 up to 50 
Mj. All models were followed from the gravitational contraction 
to an approximate age of 20 Gyr. 

2.2. Non-stop calculations: from pre main-sequence to white 
dwarf models 

We used the same version of the MESA code to compute a non- 
stop evolution from the PMS to WD cooling sequences. The 
adopted chemical composition and the mixing-length parameter 
omlt are the same as those adopted for the planet models. 



2.3. Neutron star models 

The neutron star models were computed using two different 
codes. The first is a modified version of the NSCool/TOV sub- 
routines (Page & Reddy 2006; Page, Geppert & Weber 2006). 
We computed neutron star models adopting the non-relativistic 
EOS (equation of state) by Akmal, Pandharipande & Ravenhall 
(1998) (APR). 




log age (yrs) 

Fig. 1. Time evolution of T(M, EOS) for some planets with masses 
between 0.1 and 50 My. 



Another set of subroutines was used to explore the role of the 
EOS in the invariance of the function F(M, EOS). We considered 
seven different EOS that were calculated with the model from 
Hempel & Schaffner-Bielich 2010, using the following relativis- 
tic mean-field nucleon interactions: DDI (Typel et al. 2010), 
FS Ugold (Todd-Rutel & Piekarewicz 2005), NL3 (Lalazissis, 
Ktinig & Ring 1997), SFHo, and SFHx (Steiner et al. 2012), 
TM\ (Sugahara & Toki 1994), and TMA (Toki et al. 1995). 
Some subroutines were added to compute the apsidal-motion 
constants, the moment of inertia, and the gravitational potential 
energy. 

As we have shown in Claret (2012), the effects of general 
relativity on the moment of inertia and gravitational potential en- 
ergy are not important for PMS, main-sequence (MS) stars, WD 
and, as we will see below, for planets. However, for consistency, 
we used the relativistic formalism throughout. The moment of 
inertia Iqr and the gravitational potential energy Q.qr were com- 
puted using the following equations (for the approximation of 
the moment of inertia, see Ravenhall & Pethick 1994 and for the 
gravitational potential energy, see Misner, Thorne & Wheeler 
1973): 



Jgr - 
Igr — 



y J A(r)r 4 [p(r) + P(r)/c 2 ]dr, 



J0R =(fi G RR) 2 M 



V + «v ) 

Q cs = -4tt JT r 2 p [A 1/2 (r) -\\dr = -a CR ^ 



(1) 

(2) 



where P(r) is the pressure, p(r) the energy density and the aux- 
iliary function A(r) is given 

by [l - ^Tpr 1 ] ■ Equations (1) and 
(2) were integrated through a fourth-order Runge-Kutta method. 
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Fig. 2. First panel: time evolution of T(M, EOS) for some stellar models evolving from the PMS to the WD stage. The mass range is 0. 1-1 .7 M Q 
(see symbols in the left upper corner). The adopted initial chemical composition is X = 0.73, Z= 0.01 and the mixing-length parameter a MLT = 
1.5. Second panel: specific nuclear power for the same models. 



3. Gaseous planets 



Figure 1 displays the calculations of F(M, EOS) for gaseous 
planets with masses ranging from 0.1 up to 50 Mj. Although 
A(R) tends to 1 .0 for all gaseous planet models, for consistency, 
we kept the definition of F(M, EOS). The refinement of the expo- 
nent in the above function is explained in Sect. 5. The evolution 
of each model was followed from the gravitational contraction 
to an approximate age of 20 Gyr. 

An inspection of Fig. 1 shows that Y(M, EOS) is constant 
for all evolutionary phases (« 0.4). This is the same value 
we obtained for PMS, WD and NS models (Claret 2012). We 
performed some tests by varying the initial chemical compo- 
sition and omlt and the results are independent of the input 
physics. With these calculations we extended the invariance of 
T(M, EOS) to another category of celestial bodies, the planets. 

As we show below, the invariance of the aforementioned 
function during all planetary evolutionary phases is connected 
with the specific nuclear power that, for the gaseous planets we 
analysed here, is below a threshold that is determined in the next 
section. 



4. Connection between the specific nuclear power 
e N and T(M, EOS) for PMS to the WD cooling 
sequences 

In a previous paper (Claret 2012) we argued that the drastic 
changes in a/3 during the MS, thermally pulsing asymptotic giant 
branch, and "blue loops" for PMS-WD models could be related 
to the presence/absence of nuclear reactions, with chemical in- 
homogeneities, etc. To try to elucidate this point we computed 
some PMS-WD models with masses varying from 0.1 up to 1.7 
M . We found that the specific nuclear power ejv is directly cor- 
related with the variation of Y(M, EOS), as shown in Fig. 2. The 
first peaks in for the models with 1.7 and 1.1 M just before 
the ZAMS (zero age main sequence) are due to the reduction 
in the original I2 C content through the nuclear reactions I2 C(p, 
y) 13 N(/3 + , v) 13 C(p, y) 14 N. During the PMS stage the function 
is conserved for all models and also at the final stages of WD, 
regardless of the initial mass of the model. Because the PMS 
duration depends on the initial masses, F(M, EOS) begins to in- 
crease for the most massive models first. For the least massive 
model the function is conserved almost throughout, except dur- 
ing a very short interval (log age « 12.4), but the PMS value is 
recovered a little later. 

By inspecting Fig. 2 we can deduce that there is a thresh- 
old below which F(M, EOS) is invariant. This limit is around 
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NL3 

SFHo 

SFHx 

TM1 
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Table 1. Neutron stars and EOS: exponents for the function A(R) 
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Fig. 3. Mass-radius relation for neutron stars. The adopted EOS 
are identified by the symbols displayed at the right upper corner of the 
figure. Only stable models are shown. 



10~ 3 ergs s _1 g -1 (indicated by horizontal dashed line in the sec- 
ond panel of Fig. 2). Above this threshold, the stellar mass 
distribution changes so drastically due to the presence of nu- 
clear reactions and the subsequent chemical inhomogeneities 
that V(M, EOS) increases up to four orders of magnitude (not 
shown in the figure for the sake of clarity). This occurs mainly 
during some post-main sequence phases. The behaviour of the 
configurations during these phases can be understood with the 
help of the virial theorem. If the accelerative changes in the mo- 
ment of inertia are not large (|4 ~ 0), we have 2 E - Q = 0, 
where E is the total energy and Q is the generic gravitational 
potential energy given by aGM 2 /R. The form-factor a gives 
a measure of mass concentration of a model, in a similar way 
as the apsidal-motion constant does. As the core contracts, a 
increases and the external layers must expand to keep the gravi- 
tational potential energy constant. On the other hand, below the 
threshold the configurations are changed more smoothly and this 
function remains constant (« 0.4) as in the cases of PMS, WD, 
NS, and for the gaseous planets studied in the previous section. 
Our criterion for quasi-static or stable evolution visible in Fig. 2 
differs substantially from that used by Quiroga & Claret (1992a) 
(see their figures 1, 3, and 4), who adopted 0.38 < a/3 < 0.44 
from the analytical approximations by Ferronsky et al. (1978). 
In addition, the criterion by Quiroga & Claret (1992a) does not 
relate the time evolution of a/3 with the nuclear reactions or with 
the chemical inhomogeneities. 



5. Gravitational potential energy and the moment of 
inertia of neutron stars: the role of the equation 
of state and the macroscopic stability criterion 

The invariance of [ayS] GR /A 9 (R) for NS was shown in (Claret 
2012) but using only two EOS (Akmal, Pandharipande & Raven- 
hall 1998, Shen et al. 2011). Here, we extended these calcula- 
tions to seven other EOS formulations (see subsection 2.3 and 
Fig. 3) to test whether this function is independent of the EOS 



EOS 


M max ( M e ) 


exponent 


X z /N 


DDI 


2.4223 


0.79 


6.59E-6 


FS Ugold 


1.7392 


0.86 


2.13E-5 


NL3 


2.7911 


0.78 


2.76E-5 


SFHo 


2.0587 


0.82 


1.25E-5 


SFHx 


2.1301 


0.79 


2.66E-5 


TM1 


2.2130 


0.82 


4.18E-5 


TMA 


2.0217 


0.82 


6.52E-5 


All EOS 




0.80 


4.31E-5 



and also to refine the exponent of this function. First, we anal- 
ysed each EOS separately and searched for the best exponent to 
the function A(R) by means of the least-squares method. Next, 
we considered all EOS as a whole and repeated the same proce- 
dure. The results are listed in Table 1 where the merit func- 
tion x 2 IN is also tabulated, with N as the number of points. 
All resulting exponents are close to 0.80, the strongest devia- 
tion comes from FS Ugold, which predicts 0.86. However, this 
EOS also predicts an M max of 1 .74 M@, which contradicts with 
the value of the highest inferred NS mass (1.97 M Q ) obtained 
by Demorest et al. (2010). Probably this EOS should be ruled 
out, but we decided to keep it in our calculation of the average 
exponent (see the last line of Table 1). 

The individual values of the exponent of A(R) as well as the 
average one obtained taking into account all EOS together (Table 
1) allow us to consider that the function F(M, EOS) is invariant 
for NS and independent of the adopted EOS and of the stellar 
mass. For WD and NS as the final products of stellar evolution, 
we have shown that they recover the initial value of T{M, EOS) 
characteristic of the PMS (« 0.4), i. e., stars begin and end then- 
lives in different ways but they recover this fossil function at the 
last stages of stellar evolution. Whether this remains true for 
black holes has to be examined in future studies. 

As we have seen, the function F(M, EOS) is conserved for all 
NS models/EOS but [ffySJcp. is no longer constant and is a func- 
tion of mass (Fig. 4). Some interesting characteristics can be 
inferred from this figure. First, it can be noticed that for very low 
NS masses the effects of general relativity are not important and 
[ayS] GR tends to the Newtonian value. The second characteristic 
is related to the stability of NS models. All models displayed in 
Fig. 4 are in equilibrium, but this does not necessarily mean that 
they are stable. The main point concerning equilibrium configu- 
rations is that they may correspond to a maximum or minimum 
of energy with respect to compression/expansion. For NS the 
necessary (but not sufficient) condition of stability is given by 
(e. g., Harrison et al. 1965) 



dM 

dp c 



>0. 



(3) 



An equivalent criterion can be derived considering the clas- 
sical order-of-magnitude relations for the pressure and density, 
P a GM 1 /R A and p « M/R*, respectively. From these relations, 
we obtain (Cox & Giuli 1968) 



M 



i 3/2 



Gp 4 / 3 / 



2 (d\\\M\ 



d\nP\ 
d\np J s 



3 3 



(4) 
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Fig. 4. MsLSS-[af3] GR relation for neu- 
tron stars, including unstable mod- 
els. The two arrows indicate the 
direction of increasing central den- 
sity and the full squares denote the 
M„ mx . The same symbols as in Fig. 
3. Lower right corner inset: mass- 
energy distribution for a stable (con- 
tinuous line) and unstable (dashed 
line) model with similar masses us- 
ing DDI EOS. 




1.2 1.4 1.6 U 
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Fig. 5. Mass-a^ relation for neutron stars. The same symbols as in 
Fig. 3. The horizontal full line indicates the average value of [o^GRmax- 



where the index S refers to the entropy and Fi is an average value 
over the whole model. From Eq. 4 we derive the condition that 
a model will be stable for F; > 4/3. 

Now we use the Jacobi virial theorem to check the stability 
and derive the dimension of a NS without the need to use numer- 
ical calculations like the integration of the TOV equations. We 
mainly follow the work by Fowler (1966). It can be shown that 
if we use the Jacobi stability criterion, > 0, we obtain the 
inequality 



R 



/(4£ ri-1+26 5-3rx)< 

Km v ' 



( 3 3TV-4 
2(5 - n) 



where R$ = 2GM/c is the Schwarzschild radius, is the 
radius of a neutron star, the factor 3/(5 - n) comes from the 
gravitational potential energy for a polytrope of index n, and 
f„ » 5.07/(5 - n) 2 . This function is connected with the post- 
Newtonian expansion of the proper internal energy and the grav- 
itational potential energy of polytropes. Using the results of Eq. 
(4), FT = 3/2, ft = 0.31, & = 0.55, and n = 2 we obtain 
Rn > 4-JRs- Introducing numerical values for 1 M Q , we ob- 
tain that a typical NS must have a radius of the order of 14 km, 
which agrees well with the average radius inferred from Fig. 3. 
Indeed, for the more elaborated calculations shown in Fig. 3, the 
radius of a 1 M e NS model is between 12 and 15 km. 

Next, we derive heuristically an alternative macroscopic cri- 
terion of stability for NS based on the properties of their gravita- 
tional potential energies and moments of inertia through [ffySJcp.. 
We refer to Fig. 4, where several NS models are displayed for 
seven EOS, including unstable models. In this figure the arrows 
indicate the direction of increasing of p c and the limiting masses 
are denoted by full squares. We first consider the configurations 
to the left of the turning points. If a given configuration in this 
region is compressed (higher p c ), the resulting configuration will 
present a deficit of gravitational mass if compared with its orig- 
inal position. The gravitational force will act to return it to the 
original position. Mutatis mutandis, if we decrease p c , the star 
also returns to the initial position. For models to the right of the 
turning points, this force will tend to move them from the origi- 
nal position. Therefore, we introduce a macroscopic criterion of 
stability: the models will be stable only for positive slopes 



dM 



> 0. 



(6) 



GR 



(5) 



To illustrate this criterion of stability we plotted at the lower 
right corner of Fig. 4 the mass-energy distribution for two NS 
models with similar masses (2.23 M ) but located to the left (sta- 
ble) and to the right (unstable) of the turning point. The adopted 
EOS was DD2. The steeper profile near the centre correspond- 
ing to the model located to the right of the turning point indi- 
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cates that a perturbation in the star probably induces it to col- 
lapse gravitationally to a black hole. 

As we have seen (Eq. 4), the critical value of F[ is 4/3 for 
Newtonian configurations. However, the effects of the general 
relativity on r c) , can be significant for NS, for example. Shapiro 
& Teukolsky (1983) have shown that these effects can be ex- 
pressed as 



4 GM 



Rc 2 



(7) 



where a is a positive constant of the order of unity. Therefore, 
the effect of the General Relativity is to increase the value of 
T cr i and making NS more unstable. On the other hand, we have 
shown in the beginning of this section that [o)6]gr /A & (R) « 0.4 
for all NS models. Combining Eq. 7 with the definition of the 
function A(R), we can show that 



7T- 4 , ([a/3] 



iCR 



0.4 



[a/3] 



CR 



(8) 



Equation 8 presents some interesting aspects. The first one is 
connected with the Newtonian configurations. For PMS and 
WD sequences we have shown (Claret 2012) that [ay8] GR — > 
[a/3] 

Newt ~ 0-4. Therefore, for these configurations the relativis- 
tic corrections are very small. The other point is that which links 
Fig. 4, Eq. 6, and Eq. 8. The behaviour of [o^Jgr for NS is 
a key for understanding the stability of these stars and also for 
determining M max . 

The general features shown in Fig. 4 as well as the alter- 
native criterion of stability are due to the high compactness of 
NS and are directly related to General Relativity. As mentioned 
before, new observational data on NS masses gradually act as a 
discriminator and some EOS should be discarded. At our present 
level of knowledge, the EOS FS Ugold must be ruled out. In the 
framework of Fig. 4 the limiting masses M max are located in 
a narrow range (0.72-0.85) of [qt6] gr . In Fig. 5 we show the 
behaviour of »gr as a function of mass. By inspecting this fig- 
ure we note that ccgr is almost independent of the EOS at the 
turning points, i.e., for the M max . The value of aGRmax obtained 
with a least-square fitting is 1.15: the strongest deviation is 0.05. 
Surely, acRmax can be refined with the contribution of improved 
EOS and new observational data of NS. 

In the present investigation, we only included relativistic 
EOS models and, in part the results for the non-relativistic APR 
EOS. The seven relativistic mean-field EOS of Hempel et al. 
cover a broad range of possible neutron star mass-radius rela- 
tions. Nevertheless, it would be good to include additional non- 
relativistic EOS and maybe also EOS with exotic degrees of free- 
dom such as hyperons or quarks, to further validate the model - 
independency of the function T(M, EOS). 



6. Summary 

We have investigated some aspects of the stellar evolution con- 
cerning the Jacobi virial equation, the function F(M, EOS), and 
the stability criterion. We summarise here the main results: 
1) We analysed the internal structure of gaseous planets with 
masses varying from 0. 1 up to 50 Mj. We showed that the func- 
tion F(M, EOS) is invariant for all models throughout the plan- 
etary evolution. In this way we extended the invariance of the 
mentioned function to gaseous planets. 



2) We found a connection between the strong variations of 
T(M, EOS) during the intermediary evolutionary phases with the 
specific nuclear power. We also found a specific nuclear power 
threshold. Below this limit the function is invariant (« 0.4) for 
PMS -WD models. 

3) It was shown that the function F(M, EOS) for NS models is 
independent of the EOS and of the stellar mass. 

4) A macroscopic stability criterion for NS models was intro- 
duced on the basis of the relativistic product [q^Jgr- 

5) Although the observational data for NS does not allow us to 
distinguish among all EOS yet, we determined that it seems to 
be a single maximum value of [a] GR for the turning points. 

6) We confirmed that regardless of the final products of the 
stellar evolution, NS or WD, they recover the initial value of 
[a/3] GR /A % (R) ~ 0.4 acquired at the PMS. The case of black 
holes will be subject of a future study. 

Acknowledgements. We thank the anonymous referee for his/her useful com- 
ments and suggestions. We would like to thank B. Rufino and V. Costa for their 
comments. The Spanish MEC (AYA2009-10394, AYA2009-14000-C03-01) is 
gratefully acknowledged for its support during the development of this work. 
M. H. acknowledges support from the High Performance and High Productivity 
Computing (HP2C) project, and the Swiss National Science Foundation (SNF) 
under project number no. 200020-132816/1. M. H. is also grateful for sup- 
port from ENSAR/THEXO and CompStar, a research networking program of 
the ESF. This research has made use of the SIMBAD database, operated at the 
CDS, Strasbourg, France, and of NASA's Astrophysics Data System Abstract 
Service. 



References 

Akmal, A., Pandharipande, V. R., & Ravenhall, D. G. 1998, Phys. Rev. C, 58, 
1804 

Claret, A. 2006, A&A, 453, 769 

Claret, A. 2007, A&A, 467, 1389 

Claret, A. 2012, A&A, 543, A67 

Claret, A., & Gimenez, A. 1989, A&AS, 81, 37 

Demorest, P. B., Pennucci, T., Ransom, S. M., Roberts, M. S. E., & Hessels, J. 

W. T. 2010, Nature, 667, 1081 
Cox, J. P., & Giuli, R. T. 1968, Principles of Stellar Structure, Gordon & Breach, 

2 

Ferronskii, V. I., Denisik, S. A., & Ferronskii, S. V. 1978, Celest. Mech., 18, 113 
Fowler, W. A. 1966, ApJ, 144, 180 

Harrison, B. K., Thorne, K. S., Wakano, M., & Wheeler, J. A. 1965, Gravitational 

Theory and Gravitational Collapse (University of Chicago Press) 
Hempel, M., & Schaffner-Bielich, J. 2010, Nucl. Phys. A, 837, 210 
Lalazissis, G. A., Konig, J., & Ring, P. 1997, Phys. Rev. C, 55, 540 
Misner, C. W., Thorne, K. S., & Wheeler, J. A. 1973, Gravitation, ed. W. H. 

Freeman (San Francisco) 
Page, D., & Reddy, S. 2006, Annu. Rev. Nucl. & Part. Sci., 56, 327 
Page, D., Geppert, U., & Weber, F. 2006, Nucl. Phys. A, 777, 497 
Paxton, B., Bildsten, L., Dotter, A., et al. 2011, ApJS, 192, 3 
Quiroga, R. J., & Cerqueira, A. H. 1992, EMP, 58, 215 
Quiroga, R. J., & Claret, A. 1992a, ApSS, 193, 185 
Quiroga, R. J., & Claret, A. 1992b, ApSS, 193, 235 
Quiroga, R. J., & Mello, J. C. 1992a, ApSS, 193, 247 
Ravenhall, D. G., & Pethick, C. J. 1994, ApJ, 424, 846 

Shapiro, S. L., & Teukolsky, S. A. 1983, Black Holes, White Dwarfs, and Neu- 
tron Stars. The Physics of Compact Objects (New York: Wiley) 
Shen, H, Toki, H, Oyamatsu, K., & Sumiyoshi, K. 201 1, ApJS, 197, 20 
Steiner, A. W., Hempel, M., & Fischer, T. 2012, [arXiv: 1287 . 2 184] 
Sugahara, Y, & Toki, H. 1994, Nucl. Phys. A, 579, 557 
Todd-Rutel, B. G., & Piekarewicz, J. 2005, Phys. Rev. L, 95, 122501 
Toki, H., Hirata, D., Sugahara, Y, Sumiyoshi, K., & Tanihata, I. 1995, Nucl. 
Phys. A, 588, 357 

Typel, S., Ropke, G., Klahn, T., Blaschke, D., & Wolter, H. H. 2010, Phys. Rev. 
C, 81, 5803 



Article number, page 6 of [6] 



